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ABSTRACT 

Fractal concepts have been introduced in the accretion disc as a new feature. Due to the fractal 
nature of the flow, its continuity condition undergoes modifications. The conserved stationary 
fractal flow admits only saddle points and centre-type points in its phase portrait. Completely 
analytical solutions of the equilibrium point conditions indicate that the fractal properties 
enable the flow to behave like an effective continuum of lesser density, and facilitates the 
generation of transonicity. However, strongly fractal flows inhibit multitransonicity from de- 
veloping. The mass accretion rate exhibits a fractal scaling behaviour, and the entire fractal 
accretion disc is stable under linearised dynamic perturbations. 
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1 INTRODUCTION 

Accretion processes are understood to be the non-self-gravitating flow of a compressible astrophysical fluid, driven by the external gravita- 
tional field of a massive astrophysical object, like an ordinary star or a compact object or, what is perhaps of greatest interest in the context of 
accretion studies, a black hole. The matter falling into the potential well of an accretor is commonly modelled in terms of a fluid continuum, 
and standard fluid dynamical equations, c ustomized in the N ewtonian framework of space and time, suffice well to provide a satisfactory 
description of the entire accretion process (Frank et al. 2002). While this be so, in recent times fractal modelling of astrophysical accretion 
is generating some interest (Roy 2007.: .Roy & Ray 2007.) . One astrophysical material that ought readily to lend itself to a fractal approach 
— as opposed to a continuum approach — is the interstellar medium (ISM). It is recognised that the ISM is not entirely to be seen as a fluid 
continuum, and for ma ny purposes the ISM is bel ieved to p ossess a self-similar hierarchical structure over several orders of magnitude in 
scale iLarson 1981 : Falgarone et al.lll992 : Heithausen et al.lll99&). Direct H I absorption observations and i nterstellar scintilla tion measure- 



ments suggest that the structure extend s down to a scale of 10 au ( Crovisier et alj[l983 : Langer et alJll995 ; Faison et"^ll998h and possibly 
even to sub-au scales ll Hill et aljEooi) . Numerous theories have attempted to explain the origin, evolution and mass distribution o f these 
clouds and it has been established, from both observations felmegreen & Falgarone 1996h and numerical simulations (Burkert et al. 1997 ; 
Klessen et al.lll998l ; ISemelin & CombelboOOh . that the interstellar medium has a clumpy hierarchical self-similar structure with a fractal di- 
mension in three-dimensional space. The main reason for this is still not properly understood, but it can be the consequence of an underlying 
fractal geometry that may arise due to turbulent processes in the medium. 

A comprehensive mathematical description of an accreting system — either a fluid continuum or a fractal structure — will necessarily 
have to fall within the broad domain of nonlinear dynamics. This is the principal basis of the analytical methods adopted for this study. In 
the conte xt of sphe rically symmetric accretion, a previous work furnished a self-consistent description of the hydrodynamics in a fractal 
medium l Roy||2007 ). A similar description has been provided in this work for the case of a c onserved rotational flow, which, from a general 
fluid dynamical perspective, has become quite thoroughly understood ( Chandrasekhar ' 1 98 iV In the context of accretion studies particularly, 
on many occasions such flows are devised to be an inviscid sub-Keplerian compressible flow, and as such, this model has become well- 
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estabhshed in acc r etion-related literature bv now jAbramowicz & Zurek 
Chakrabarti"l990': Kafatos & Yang" 1994'; Vang & Kafatos"l995'; 'Pariev"l99( 
2002 : Das et_al. 2003: Ray 2003: Barai et al, 2004: Das 2004: Abraham et al. 2006: Das_et al. 2007: Chaudhury et al. 2006: G oswami etal 
2007 ). As a first step, the physical processes in a fractal medium have been analysed by fractional integration and differentiation jzaslavskv 
20021 and references therein). To do so, the fractal medium has had to be replaced by a continuous medium and the integrals on the network 
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of the fractal medium has had to be appro ximated by fractio nal integrals jRen et alj|2003h . The interpretation of fractional integration is 
connected with fractional mass dimension jMandelbrolll 19831) . Fractional integrals can be considered as integrals over fractional dimension 
space within a numerical factor ( lTarasovll2004f) . This numerical factor has been chosen suitably to get the right dimension of any physical 
parameter and to derive the standard expression in the continuum limit. The isotropic and homogeneous nature of dimensionality, in the 
context of an axisymmetric flow, has also been incorporated properly. 

Having established the hydrodynamic equivalen ce, a standard mathematical path has been followed to set up the steady fractal flow like 
a simple first-order autonomous dynamical system ( Strogatz 1994 ; Jordan & SmithI 1999i) and investigate the resulting critical properties, 
especially the transonic character of the flow. A general prediction that can be made about the conserved steady fractal disc flow is that 
the equilibrium conditions will only permit the existence of saddle points and centre-type points in the phase portrait of the flow. Fully 
analytical methods have been employed to identify the exact critical point coordinates, under a pseudo-Newtonian prescription to study the 
infall process. Two very significant effects of fractal properties can be discerned in consequence of this. The first is that fractal accretion flows 
assume the properties of an equivalent continuu m, albeit one w ith a more diluted effective local flow density, a result that is also known to be 
true for the case of spherically symmetric flows ( Roy & Ravl2007ii'I The se cond effect is that fractal proper t ies have an inhib itory influence on 
the multitransonic character that is commonly associated with these flows Jchakrabartill 19891 , 1 199(1 1 1996l ; lDasll2002i |2004|) . Taken together, 
the overall effect of strongly f ractal features in a disc flow is to lend it the appearance of the spherically symmetric monotransonic flow — 
the paradigmatic Bondi 1 1952h solution. 

In the Newtonian limit the mass accretion rate in a fractal flow exhibits a scaling behaviour that is governed by the fractal dimension 
of the flow. This scaling behaviour may have many interesting astrophysical implications. One particular case, that of accretion processes 
onto a pre-main-sequence (PMS) star, should have a close connection with this fractal scaling relation. It has also been shown that the sta- 
bility of fractal flows under linearised time-dependent perturbations, is governed by the same equation pertaining to a standard perturbation 



bility or iractai flows unaer linearised time-aependent perturbations, is governed by the same equation pertaimng to a 
schem e that has been used earlier for studying continuum sub-Keplerian disc flows i Ravlboojl : IChaudhurv et al.ll2006l : 



Rav & Bhattacharjee 



20071) . Given that the boundary conditions do not vary between the continuum case and the fractal case, all steady global solutions (tran- 



sonic or otherwise) are, therefore, stable. Since fractal disc flows tend to acquire the character of spherically symmetric flows, it has be- 
come possibl^tojnvokejlwai^ and dynamic evolution of the spherically symmetric flow towards its transonic 
state <Rav & Bhattacharie3l2002l : Irov & Ravlboovb . and suggest that in a fractal disc, the drive towards transonicity would also occur under 
identical guiding principles. 

It would be important to note that the entire analytical treatment has been carried out in the pseudo-Newtonian formalism, with the 
mathematical results expressed in terms of a general potential, whose resulting force field would drive the accretion process. So most of 
the conclusions reached here will remain indepen dent of the choice of a po tential. However, the specific numerical plots which have been 
presented here, have all been derived by using the IPaczviiski & WiitalHisol) potential for a polytropic gas flow. In spite of this, the general 
claims made on the basis of these particular numerical results will s uffer no qualitative alteration, under the choic e of any other pseudo- 
Newtonian potential i Nowak & Wagonelll991 ; Artemova et alJl993) . or for an isothermal gas flow i Das et al.]l2003 ). 



2 THE HYDRODYNAMIC EQUATIONS FOR AXISYMMETRIC FRACTAL ACCRETION 

Considering a thin accretion disc that has density fluctuations over a range of length scales, it should be possible to propose that these density 
fluctuations or "dumpiness" can be approximated as a fractal structure along the radial direction in the disc. This will be especially true in the 
thin-disc regime, where there is not much scope for significant variations to occur i n the vertical dir ection IFrank et al ]|2002) . The integrals 



on this network of fractals can consequently be approximated by fractional integrals ( Ren et alj2003h , and the interpretation of this fractional 



integration can be connected to the fraction al ma ss dimension. This approximation has already been used to describe the fractal structure of 



spherically symmetrical accretion by Rovl ( 120071) . Likewise, a similar transformation of the infinitesimal length of the "clumpy" accretion 
disc is being employed here to describe it in terms of an equivalent "fractional continuous" disc. The fractional infinitesimal length for such 
a structure will be given by 

dr=flj dr, (1) 

with Zc being the inner length scale of the structure, representing the physical scale of length at which the density is correlated. The fractal 
dimension of the accretion disc is given as 2A. Hence, the mass enclosed in a thin disc of density, p, radius, r, and thickness, H, will be 

Md ^ / 27rpr dr dz ~ pffr^"^ , (2) 

Jo J-H/2 

which, quite obviously, will give the standard expression in the limit A — > 1. 

S ince there is no flow in the z direction (perpendicular to the plane of the disc), the condition of hydrostatic equilibrium i Frank et sl\ 
20021) . by balancing the pressure gradient and the gravitational force acting on a fluid element in the r~(p-z space, will give 



^ In the context of fract als, it would be interestin g to note a similar equivalence between the discrete model of diffusion-limited aggregation and the continuum 
dendritic growth model <Witten&SandeJl98lh . 
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^ (Pr Ar Acj)) Az ~ c&'p ( — \ r Ar A(j>Az 
dz \r J 



(3) 



in which P is the gas pressure, and <I> is the gravitational potential of the central accretor that drives the flow (with the prime denoting the 
spatial derivative of the potential). In the case of stellar accretion, the flow is driven by the Newtonian potential, $(r) = —GM/r. On the 
other hand, frequently in studies of black hole accretion, it becomes con venient to use an effective pseudo-Newtonian potential that imitates 
gener al relativistic effects in the Newtonian construct of space and time i Paczviiski & Wiitalll98d : Nowak & Wagone3ll991 ; Artemova et al. 

Along with the definition of a polytropic equation of state i Chandrasekhai 1939h . P — Kp^ (with the polytropic expone nt, 7, having 
the range 1 7 ^ 5/3), equation ^ can be simplified to give the condition for the thin-disc approximation Jprank et alJlOoA as 



H ■ 



T 
7$' 



1/2 



(4) 



in which the local speed of sound, Cs, has been defined by Cg = dP/dp. This definition, as well as the form of H in equation l|4}, will be 
instrumental in determining the mass density balance (the continuity condition) in a conserved fractal accretion disc. For an infinitesimal 
volume of the fractal medium, AV, this balance will be given by 



d 



dr 



AV 



Ar 



Making note of the connection that AV = Hr Ar ( 
derived as 



as 1 

dt ^ r 



d 



0, 



(5) 

for the axisymmetric thin disc flow, a more informative form of equation {5]( can be 



(6) 



2A-1 Qj. 

in which E is the surface density, defined by E ~ pH jFrank et al ] |2002h . Using the thin-disc approximation, as given by equation l|4), one 
can recast equation ^ further as 



d_ 

dt 



[^(7+l)/2j 



+ 



/W_d_ 



(7 + l)/2 T 



0, 



(7) 



where a = 2 A — (1/2). The foregoing expression gives the mass density balance equation (the continuity equation) for the axisymmetric 
fractal flow. 

Similarly, in the infinitesimal volume, AV, the balance of momentum density will imply 
^ (pvdF) =Fg+Fcf + Fp, (8) 

with V being the velocity vector, Fg and Fcf being, respectively, the gravitational and centrifugal forces acting on the mass contained in the 
infinitesimal volume, AV, and Fp being the total surface force due to the pressure acting on the full surface bounding the volume, AV. The 
first two forces have radial components only, and their magnitudes are given by 



Fg = -^'pAV 
and 

\2 _ 

Fcf = ^pAV, 



(9) 



(10) 



respectively, with A in the latter being the constant specific angular momentum of the conserved disc flow 1 Chakrabarti 1989I . 199(1 19961) . 
On the other hand, Fp has r, (ji and z components. However, only the radial component is relevant for the axisymmetric thin- disc flow. 
Any v olume element, AV, experiences a force, — dV(VP), due to the pressure exerted on it by the surrounding medium (Landau & Lifshitz 
ll987D . Translating this effect onto an infinitesimal volume element of the fractal medium, the magnitude of the radial component of the force, 
Fp, can be set down as 

Fp = TTP'^^- 

p or 

Now going back to the left-hand side of equation {§](, the total change of radial momentum is extracted as 



At 



(pvdl/) 



dv dv 
dt dr 



pAV, 



(11) 



(12) 



a result that could be derived by invoking the continuity condition, as it is given by equation And so, making use of equations {9}, JlOt . 
dllb and J12b in equation {8}, it becomes possible to obtain the final momentum balance condition for the conserved fractal disc flow as 



dv dv 1 dP ^ , , > 
dt dr p dr 



0, 



(13) 
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which actually gives a local conservation law and, as expected, has exactly the same form as that of the continuous medium. The fractal flow 
can now be described completely by equations (O and l ll3t . both of which describe the dynamics of the fields v(r, t) and p(r, t), with the 
pressure, P, having been defined as a function of p. 



3 THE STATIONARY FRACTAL FLOW AND ITS FIXED POINTS 

It is a standard practice to consider a steady flow while studying accreting systems i Frank et"al]|2002h . The two corresponding stationary 
equations which determine the drift in the radial direction, can be obtained from equations ([7} and J13t as 



(7+i)/2 r 



= 



(14) 



(15) 



_d_ 

dr L' 

and 

dv IdP , 
V— + -— +$'(r) 
dr p dr i - 

respectively. The pressure, P, is prescribed by an equation of state for the flow ( Chandrasekhar 1939h . As a general polytropic it is, as 
usual, given as P = Kp'', while for an isothermal flow the pressure is given by P = pKT/ pmn, in all of which, if is a measure of the 
entropy in the flow, 7 is the polytropic exponent, k is Boltzmann's constant, T is the constant temperature, mn is the mass of a hydrogen 
atom and p is t he reduced m ass, respectively. Since the local speed of sound is defined as Cg = {dP/dpY^^ , the way in which P has 
been prescribed jFrank et al ][2002) will affect the transonic features of the flow, because transonicity, and all its related critical aspects, are 
measured by scaling the bulk flow velocity with respect to the local speed of sound. So the exact form of P has distinctive consequences, and 
in what follows, the flow properties will be taken up separately for the two cases, i.e. polytropic and isothermal. 



3.1 Polytropic flows 



With the polytropic relation specified for P, it is a straightforward exercise to set down in terms of the speed of sound, Cs, a first integral of 
equation dlSb as. 



V A 
^ + nc,2+$(r) + — 



(16) 



in which n = (7 — 1) ^ and the integration constant £ is the Bernoulli constant. The first integral of equation J14l ( could similarly be 
obtained as 



(17) 



where M = {'yK)"rh Jchakrabartilll99dl 19961) with rh, an integration constant itself, being physically the matter flow rate. 

To obtain the critical points of the flow, it should be necessary first to differentiate both equations l |16t and l |17t . and then, on combining 
the two resulting expressions, to arrive at 



^2 2x d 2n 2ii^ 

/ dr r 



r$' + -13'^cl [2a -r 



(18) 



with 0^ =2(7 + 1)^^. The critical points of the flow will be given by the condition that the entire right-hand side of equation dlSt will 
vanish along with the coefficient of d{rP')/dr. Explicitly written down, following some rearrangement of terms, this will give the two critical 
point conditions as 



rc*'(rc) 



2cr 



$"(rc) 



$'(rc) 



(19) 



with the subscript "c" labelling critical point values. 

To fix the critical point coordinates, Vc and Tc, in terms of the system constants, one would have to make use of the conditions given by 
equations l |19t . along with equation dl6t . to obtain 

-1 



27 



7- 



1 



A" 



2(7 - r< 



$"(r 



$'(rc 



2r? 



(20) 



from which it is easy to see that solutions of rc may be obtained in terms of 7, A, a (or A) and £ only, i.e. rc = fi (7, X,a,£). Alternatively, 
Tc could be fixed in terms of 7, A, a and A4. By making use of the critical point conditions in equation illi . one could write 

n 2{n + l) 



47r^/3Vg'" 
7$'(rc) 



2 



rc$'(r-c) 



2a 



$'(rc) 



(21) 
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whose obvious implication is that the dependence of will be given as rc = 72(7, A, cr, A^). 

The slope of the continuous solutions which could possibly pass through the critical points are to be obtained by applying the L' Hospital 
rule on equation (Tsj, at the critical points. This will give a quadratic equation for the slope of the stationary solutions at the critical points 
themselves in the r — phase portrait. The resulting expression will read as 



+ ^0 = 0, 



in which the constant coefficients, Z\ and Zi^, are given by 



^1 = 
and 

2o = 



7 



7 \7 + 1 



2(7 - r, 



$"(rc) 



$'(rc) 



6A^ 



+ 2<l>"(rc) + 



7+1 



2cr d 
r? dr 



+ 



7+1 



2a 



(22) 



$"(rc) 
" <E>'(rc) 



respectively. 



3.2 Isothermal flows 



For isothermal flows, a linear dependence between P and p is applied in equation J15t , as the appropriate equation of state. On doing so, the 
first integral of equation J15b is given as 

„2 y2 



y + c,^lnp+$(r) + — =C, 



(23) 



with C being a constant of integration. For flow solutions which specifically decay out to zero at very large distances, the constant C can be 
determined in terms of the "ambient conditions" as C = Inpoo- The thickness of the disc, H, is determined simply by setting 7 = 1 in 
equation l|4j, and in a likewise manner the first integral of equation will imply the continuity condition as 



2 2 2ct 

p V r 



$' 47r2c| 

As has been done for polytropic flows, both equations i23t and l l24t are to be differentiated and the results combined to give 



2N d 2^ 2^;^ 

Od;:(-) = — 



+ -ct 



from which the critical point conditions are easily identified as 

-1 



\2 

rc"I>'(rc) 5- 

ri 



2(7 - r, 



(24) 



(25) 



(26) 



$'(rc). 

In this isothermal system, the speed of sound, Cs, is a global constant, and so having arrived at the critical point conditions, it should be easy 
to see that rc and Vc have already been fixed in terms of a global constant of the system. The speed of sound can further be written in terms 
of the temperature of the system as Cs — QT^^^ , where — (K/pmn)^''^, and, therefore, it should be entirely possible to give a functional 
dependence for Tc, as rc — f-i{\, (7, T). The slope of the stationary solutions passing through the critical points in the r — v'^ phase portrait 
is obtained simply by setting 7 = 1 in equation l l22t . 



4 THE FRACTAL ACCRETION DISC AS A DYNAMICAL SYSTEM 



The equations governing the flow in an accreting system are all fluid dynamical equations , and, in the kind of i nviscid regime chosen 
for this study, fall under the general category of first-order nonlinear differential equations ( Jordan & SmithI 19991) . However, there is no 
standard prescription for deriving rigorously analytical solutions of these equations. In this situation, while a numerical integration of the 
flow equations is most often the only recou rse, an alternat i ve approach could als o be adopted by setting up the governing equations to 
form a standard first-ord er dynamical system ( istrogatJl994l : ljordan & Smithlll999l) . This is a very usual practice in general fluid dynamical 
studies I Bohr et al .Il993h . and avoiding involved numerical processes, this approach allows one to derive significant physical insight about the 
behaviour of the flows. To take a first step towards this end, for the stationary polytropic flow, as it has been given by equation (I18t, it should 
be ne cessary to parametrise this equation and set up a coupled autonomous first-order dynamical system as Jstrogatdl 1994 1 Jordan & SmithI 
I999I) 



d , 2s 

dr 
d7 



— r [V 



(27) 
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with r being an arbitrary mathematical parameter. Apropos of accretion studies, this kind of parametrisation has been reported befo re jRav & Bhattachariee 
2OO2I: kfshordi & Paczynskilliooi : Ichaudhurv et ahlbood : iRav & BhattachariedboOTl : iMandal et alj2007l:lGoswami et alJlOOTl). Some ear- 



lier works in accretion had also m ade use of the general mathematical aspects of this approach jMatsumoto et al.lll984f : lMuchotrzeb-Czemv 
1986l:lAbramowicz & Katolll989h 



The critical points have themselves been fixed in terms of the flow constants. About these fixed point values, on applying a perturbation 
prescription of the kind = vl -\- , Cg = c^c + and r = rc + 5r, one could derive a set of two autonomous first-order linear 
differential equations in the 5r — 5v^ plane, with Sc^ itself having to be first expressed in terms of 5r and 5v^, with the help of equation ( 117b 
as 



&4 

_2 



7 - 



7 + 1 



5f_ 



2(7 — rc 



$"(rc) 
$'(rc) 



(28) 



The resulting coupled set of linear equations in 5r and Sv^ will be given as 



1 d 

2t;2 dr 

j__d_ 

rc dr 

in which 



{5t) 



7- 



A 

2 V7+I 
27 



5v 



2A" 



+ *'(rc) +rc'l>"(rc 



2 $'(rc 



) 2 „ 

-CscB + — 



7+1 



5r 



7 + 1 



Sv-" -A 



7 + 1 



5r , 



(29) 



A^r, 



$"(rc) 



2<T, 



Trying solutions of the kind Sv'^ 
stability matrix, as 



<E>"'(rc) $"(rc) 



$"(rc) " $'(rc) ■ 

exp(fir) and Sr ~ exp(f2r) in equations i29h will lead to an expression for the eigenvalues, 51, of the 



^2 ^ 4rc$'(rc)cf, 



(7 + 1)^ 



{'y-l)A-2'y{2a + l + A) + 2-yB 



2a\\ 



47 + (7 - 1) ^ + 27S 



(30) 



where X^ir) = r^$'{r). 

For isothermal flows, starting from equation ([25}, a similar expression for the related eigenvalues may also be derived, and it can be 



shown that this will correspond simply to the special case of 7 = 1 in equation ( I30t . Once a critical point coordinate, (rc, v^), has been 
ascertained, it becomes an easy task to find the nature of that critical point by using rc in equation i30i . Since it has been discussed in 
Section[3]that rc is a function of 7, A, a and T for isothermal flows, and a function of 7, A, a and £ (or A4) for polytropic flows, it effectively 
implies that can be expressed as a function of the flow parameters for either kind of flow. A generic conclusion that can be drawn about the 
critical points from the form of fi^ in equation ( I30t . is that for a conserved pseudo-Schwarzschild fractal disc flow, the only admissible critical 
points will be saddle points and centre-type points. For a saddle point, fi^ > 0, while for a centre-type point, < 0. Once the behaviour of 
all the physically relevant critical points has been understood in this way, a complete qualitative picture of the flow solutions passing through 
these points (if they are saddle points), or in the neighbourhood of these points (if they are c entre-type points), can be constructed , along with 
an impression of the direction that these solutions can have in the phase portrait of the flow ('strogat zll994 ; Jordan & Smith 1999). An earlier 
study has shown how the critical point behaviour (and multitransonicity) varies with flow parameters like A and £ tChaudhurv et al...2006l) . 
Similar aspects of the accretion disc will be considered now with respect to the fractal properties in the flow. 



5 MULTITRANSONICITY IN RELATION TO FRACTAL FEATURES 



Multitransonicity in the flow is possible when there are more than one critical point through which a continuous solution will pass, connecting 
the event horizon of the accretor to the outer boundary of the flow. Saddle points are most likely to fulfil this requirement, and so multitran- 
sonicity would imply the existen ce of more than one saddle point in the phase portr ait of the flow. Disc flows, conserved or otherwise, have 
been known to be multitransonic jLiang & ThomsorJI 1 98d : lchakrabarti|[l99ol 1 19961) with a variety of implications, most notably in regard to 



shock formation (I chakrabartilll99d.ll99d IPaslbooj) . So a clear knowledge of the number of critical points in the flow, and their nature, is 



essential to understanding multitransonicity. 

For a general polytropic flow, the position of the critical points can be determined by extracting the roots from either equation ( I20t 
or ( 121b . after prescribing an explicit functi onal form fo r the pseudo-Newtonian potential, $(r). Various presc riptions are used for "l'(r), 
depending upon specific necessities (Paczvnski & Wiita 1980: Nowa k & Wag oner 1991.: .Artemova et al. 19961) . but all of these potentials 
have the common purpose of describing general relativistic effects in a semi-Newtonian framework, by avoiding the dif ficulties of a purely 
genera l relativistic treatment. One such potential that is very regularly invoked in accretion-related literature, is due to iPaczvnski & Wiita 
(liiB). t is given as 



$(r) 



2(r-l)' 

in which the radial distance, r, has been scaled by the Schwarzschild radius, Vg 



which follows the standard definition of r^ 



(31) 

2G'M/c^ with 
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M being the mass of the black hole, G the universal gravitational constant, and c the velocity of light in vacuum. All relevant velocities will 
be scaled with respect to c, and for co nvenience one may set G = c = M = 1. All derived quantities like energy and angular momentum 
will be scaled accordingly ( D asl2002 ). 

Using the functional form of given by equation i3U , in equation J20t . will lead to a quartic polynomial in the standard form 



rt + 2Arl + Brl + 2Cr^ + D = 0, 
in which 



(32) 



A = 



C = 



4£ {a + 1) 
■iS {a + 1) 



a + 1 



7 



2(7 + 1 - 



7-1 
47 



2£ (2a + 1) 



B 



7 - 



D = 



2£ (a + 1) 

A" r 27 



27 

7-1 



a{2£- 



7- 



1 



2£ (a + 1) 

It is not difficult to see that equation i32\ will yield four roots. These roots can be found completely analytically by using Ferrari's method 
for solving quartic equations. In order to do so, a term going like (arc + b)^ is to be added to both sides of equation i32l , and then 
the resulting left hand side is required to be a perfect square in the form (r^ + Arc + C)^' that the full equation will be rendered as 
(r^ + Arc + C)^ = ("'"c + b)^- This will give three conditions going as 



= D + fe^ 



2 _ {AC - cY , _ AC - c 

Eliminating a and b, will deliver an auxiliary cubic equation in C,, going as 

2^ ~ bC + 2 [AC - D)C+ {BD - A^D - C^) = , 

which, under the transformation, C = rj + (_B/6), can be reduced to the canonical form of the cubic equation, 

+ Pri + Q = 0, 
with 



(33) 



(34) 



(35) 



B' 



+ {AC - D) 



— ^j{AC-D) + ^{BD~ A^D 



12 ' ' ' ^ 108 

Completely analytical solutions for the roots of equation J35t can be obtained by the application of the Cardano-Tartaglia-del Ferro method 
for solving cubic equations. This will lead to the solution 



1/3 



+ 



1/3 



(36) 



with the discriminant, T>, having been defined by 

D = — H . 

4 27 

The sign of T) is crucial here. If D > 0, then there will be only one real root of r) given directly by equation ([36). On the other hand, if 
2? < 0, then there will be three real roots of r], all of which, under a new definition. 



(37) 



I? = arccos 



-Q/2 



can be expressed in a slightly modified form as 

+ 27r (j - 1) 



Tjj ~ 2\j COS 



(38) 



(39) 



with the label j taking the values j — 1, 2, 3, respectively, for the three distinct roots. 

Quite evidently, the sign of D and, in consequence, the number of real roots of ri will be determined by the value of the flow parameters 
like £, A, 7 and A (on which a depends). For fixed representative values of three of these parameters — £, X and 7 — the continuously 
varying dependence of ?7 on A has been shown in Fig.[T] The two distinct regions in the plot, corresponding to 7? > and 7? < 0, show 
the existence of a single and triple real roots, respectively. The one root that runs continuously across both the regions is the one that will be 
relevant for all subsequent analysis. 

Once 77 is known, it is a simple task thereafter to know the value of a and 6. With this having been accomplished, all the four roots of 
Tc from equation l l32t can be identified by solving the two distinct quadratic equations. 



Arc + C = ± {arc + b) 



(40) 



corresponding to the two signs on the right hand side. Each of these roots will determine the spatial position of a critical point in the flow 
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A 

Figure 1. For the parameter values, E = 0.0001, 7 = 4/3 and A = 1.7, the roots of?? are divided into two classes, spanning across the range, 0.5 A ^ 1. 
The segregation of the two classes is marked by the vertical line passing through the point, P, which keeps shifting to the left of the plot, as A is increased. 
To the left of this line is the region where 2? > 0, with 77 having one root only. On the right is the region where © < 0, with 77 having three real roots. The 
globally relevant root is given by the locus at the bottom, running continuously through both the regions. The cusp on the right actually comprises the loci of 
the other two roots of r?. These two distinct roots merge at P, where A ~ 0.62. 



space. Of these four roots, one is always real, and it is always located within the event horizon of the central black hole. Hence, it will be 
irrelevant for this study. Of the other three roots, one remains always real, while the other two, depending on the flow parameters, can only 
be either both real or both complex (complex roots always occur in pairs). In the former case there will, therefore, be three real roots on 
physically relevant length scales outside the event horizon, and this situation will correspond to multitransonicity. 

Taking up the case of the root that always remains real at a given physical distance in the flow region, it can be seen from Fig.|2]that as 
A decreases (i.e. as the fractal properties become more pronounced), the position of this critical point, r^, is monotonically shifted outwards. 
The nature of this critical point can be determined from the plot in Fig. [3] which shows that > for the entire relevant range of A. This can 
only mean that the critical point implied by this particular root of t-c, will always be a saddle point, and any solution passing through it, will 
necessarily be transonic. Practically speaking, all of this is just how it should be. To make accretion a feasible process, there has to be at least 
one solution that will link the outer boundary to the event horizon of the black hole. This can be achieved only if the first critical point that 
a fluid element encounters while travelling towards the accretor, after having started from an outer boundary region, is a saddle point. So the 
steadfast existence of this particular saddle point (which is the outermost critical point in the flow region), preserves the conditions necessary 
for transonicity to develop. But why does the position of this saddle point get shifted outwards, as the fractal properties become progressively 
more conspicuous? An earlier study on fractal flows in spherically symmetric geometry has addressed this very question by arguing that a 

I, and in being so, allows gravity to win over 
gas pressure resistance. The same feature has appeared in this axisymmetric case. Here transonicity occurs when gravity wins over both the 
centrifugal effects and the pressure effects in an axisymmetric flow. The conditions for transonicity are, therefore, more conducive when 
either of the two effects is weak. A relatively more dilute gas will offer a feeble pressure build-up against gravity, and so the gravitational 
pull can triumph even at greater distances. It is for this reason that when the flow is more fractal, that the transonic length scale is shifted 
outwards. The analogy with the spherically symmetric case is quite apt here. 

While this be the case for the outermost critical point, which is always a saddle point, to have an appreciation of how multitransonic 
properties get affected by the fractal nature of the flow, it shall be necessary to turn to the other two critical points — the innermost critical 
point and the middle critical point — both situated very close to the accretor. These two critical points cannot exist without each other, 
and their collective behaviour bears this out, something that is shown in Fig.|4] The positions of the innermost critical point and the middle 
critical point approach each other as A is decreased, till they coalesce. Beyond this point, multitransonicity is not possible, and the accreting 
system will only have one saddle point left in its phase portrait (i.e. it will be monotransonic). The nature of both these critical points can be 
discerned from Fig.|5] The innermost critical point is always a saddle point, and the middle critical point is always a centre-type point. The 
two critical points merge under the condition, 57^ = 0, and this merging might be viewed as the mutual annihilation of a stable centre-type 
point (with < 0) and an unstable saddle point (with > 0) in the phase space. Significantly enough, the two critical points meet when 
the lower cut-off value of A ~ 0.62. This is the same value of A, for which T>, as defined from equation i37t . vanishes, and the number 
of real roots of 77 changes from one to three. In physical terms all of this will mean that if A is reduced further, i.e. if the flow becomes 
more fractal, it will rule out the possibility of multitransonicity completely. Strongly fractal disc systems, therefore, would be devoid of the 
multitransonic character usually associated with conserved continuum disc accreting systems. 

So the general conclusion that one might draw regarding the influence of the fractal features in the flow, is that the flow is made to 



fractal medium behaves like an effectively more dilute continuous medium ( Rov & Ravll2007 
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Figure 2. The position of the outermost critical point is shifted outwards when fractal properties emerge progressively stronger (as A decreases). The parameter 
values chosen are £ = 0.0001, 7 = 4/3 and A = 1.7. 




Figure 3. For the entire relevant range of A, the outermost critical point behaves like a saddle point (Q^ > 0). The parameter values remain the same at 
£ = 0.0001, 7 = 4/3 and A = 1.7. The saddle-type behaviour of this critical point makes transonicity possible in the flow. The fractal properties of the flow 
serve to enhance this effect. 



behave like a more dilute continuum. While this helps in generating the transonic solution, in a somewhat paradoxical sense, the same fractal 
features turn out to be inimical to multitransonicity in the flow. These inferences have been made by studying a general polytropic flow, but 
arguably they will not be violated when the flow is to be considered as isothermal. In this limit also it is possible, from equation J26t . to 
derive a quartic polynomial just as in equation l l32b . with the coefficients A, B, C and D being set down, respectively, as 



4ci (a + 1) ' ci (a + 1) ' ^ cf (a + 1) ' d {a + 1) 

The rest of the mathematical analysis will remain the same, as it has been for the polytropic flow. 



6 FRACTAL SCALING FOR THE MASS ACCRETION RATE 

With the knowledge that the critical conditions in the flow can be represented entirely in terms of flow parameters like A4 (or m) and £, it 
now becomes possible to find a direct connection between the steady accretion rate and the mass of the accretor whose gravitational force 
field drives the fractal accretion flow. This should be particularly revealing for the case of the Newtonia n potential, $(r) = —GM r^^. In 



this limit it is can be shown from equation i20i that rc is scaled as GM£ ^ and A is scaled as GM£ ^''^ i Ray & Bhattacharieel2007 ). Using 
these scaling conditions in equation one can find a scaling behaviour for the accretion flow rate, m, going as 
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Figure 4. The upper arm of the cusp is the locus of the position of the middle critical point, for varying A. The lower arm is likewise traced out by the 
innermost critical point. The two tracks coalesce at Tc ~ 3.4, when A ~ 0.62. For any lower value of A, these two critical points will cease to exist, with 
there being only one saddle-type critical point, as shown in the earlier plots. The point of merging of the two critical points shifts to the left of the plot, if A is 
increased. As before, here the chosen parameter values are E = 0.0001, 7 = 4/3 and A = 1.7. 




Figure 5. The nature of the innermost critical point is indicated by the upper arm of the cusp. This critical point always behaves like a saddle point (f2 > 0). 



The middle critical point is always centre-type (Q? < 0). The two critical points merge when f!^ 
possible for a limited range if A, in this case 0.62 < A sC 1, when £ = 0.0001, 7 = 4/3 and A = 



= and A 
1.7. 



0.62. Therefore, multitransonicity is only 



(41) 



This result is particularly interesting since it indirectly provides some insight about the nature of density dumpiness or the fractal structure 
of the accretion disc in a specific astrophysical situation. 



In the case of accreting pre-main-sequence (PMS) stars, it has been argued jPadoan et al J 120051) that even if accretion takes place 
thr ough the form ati on of accretion di scs on relatively small length scales, on larger scales it can well be approximated as spherically symmet- 



Bondi 



1 1952) orlSondi & Hovle 1944 ) accretion. This scheme is in good agreement with both numerical simulations and observational 



ric 



results JPadoan et al.l 



20051). Now pre- main- sequence stars a re embedded in the interstellar medium (ISM), and numerous observational 



1 I — II 1 I II 1 I II n I II 

1985; Falgarone et al. 1992; Langer e t_alj|1995l ; lElmegreen & Falgaronelll99q ; IBurkert et alj|1997 

Faison et al. 1 9981 ; Heithausen et al. 1998. ; .Klessen et al..,1998. ; .Hill et al. .2005 ), using a variety of techniques, have directly or indirectly 
suggested that the interstellar medium has a fractal structure over a wide range of scale l engths. So, it is more realistic to consider the case 
of accretion onto PMS stars as fractal accretion. In this very context it has been shown by Roy ( 2007 ^ that the numerical simulations and the 
currently available observational res ults pertainin g to PMS star accretion rate, are consistent with the model of spherically symmetric fractal 
accretion. It has also been argued bv lRovi ( 120071) that the accretion rate is scaled as m ~ M^"\ where D{< 3) is the mass dimension of 
the surrounding ISM. So, if the proposed global model for accretion onto PMS stars — spherically symmetric accretion at large distances 
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from the star and disc accretion close to it — is correct, tiien in the steady state, going by equation l l41t . one will have 2A = D — 1 < 2. 
This will imply that the disc will also have a fractal structure, if the surrounding medium, from which the large scale accretion takes place, is 
fractal in nature. 



7 DYNAMIC ASPECTS OF THE AXISYMMETRIC FRACTAL FLOW 



The fractal nature of the accreting matter acts against multitransonicity. Nevertheless, this does not preclude transonicity itself. If anything, 
the static flow shows that prominent fractal features will certainly bring about the existence of one saddle-type critical point in the phase 
portrait of the flow. T he overall app earance of a conspicuously sub-Keplerian fractal fl ow will , there fore, be like that of a monotransonic 
spherically svmmetric iBondil J1952I) flow. This is fortunate in many respects, because the lBondil 1 119520 solution has by now become quite a 
well-understood and a regularly-invoked paradigm in accretion studies, and any close convergence between it and the fractal disc flow, will 
pave the way for making insightful comparisons. 

The most significant of these is that there will at least be one continuous solution which will connect the outer boundary of the 
flow to the event horizon of the black hole accretor, in a process that will make black hole accretion realisable in its e xpected fashion. 
The steady state conditions will imply that this solution will pass through a saddle point, which is known to be unstable ( Jordan & SmithI 
19991) . Any solution passing throu gh this point will suffer from the problem of fine tuning the outer boundary condition with infinite pre- 



cision jRav & Bhattachariedl2002h . Thes e difficulties can b e avoided through the dynamics, and indeed in the case of fractal spherically 
symmetric accretion, it has been s hown I Roy & Ravi 1 20071) that the stronger is the fractal nature of the flow, the more successful is the 
time-evolutionary drive towards the lBondil ( Il952n solution. This happens simply because a fractal medium translates equivalently as a con- 
tinuum with an effective lesser density, i.e. a fractal flow can be construed as a more dilute continuum. Since a sizeable resistance against 
gravity-driven transonicity happens due to the pressure build-up in the flow (which, through the polytropic relation, is connected to the local 
flow density), any dilution in the flowing medium will detract from the opposition against transonicity. Hence, this entire effect will enable 
an accreting solution to cross the sonic horizon smoothly, premised on the condition that this solution will also correspond to a minimum 
possible energy configu ration, and, concomitantly, to the maximum possible inflow rate i Bondil 19521 : Garlick 19791 : Rav & Bhattacharjee 



2002l : lRov&Ravll2007l) 



This line of reaso ning can be carried over full y to the case of inviscid axisymmetric accretion. As a matter of fact, in the context of 



inviscid thin-disc flows. lRav & Bhattacharjee! ( 120071) have alre ady argued that transonicity is determined and governed by the same dynamic 
and non-perturbative criteria, as they are for the Bondil ( 1952) flow. The analogy with the lBondil Il952l) solution will be more so true for 



sub-Keplerian solutions with low angular momentum (having feeble centrifugal effects pitted against gravity). Fractal features in this kind of 
disc configuration will only serve to facilitate transonicity even further. 

As opposed to a completely non-perturbative approach to the question of time-dependence in the fractal disc flow (all of which will 
require the mathematics of partial differential equations), it would also be worthwhile to study the properties of the background stationary 
flow under the influence of a linearised time-dependent perturbative effect. This will yield much information on the global stability of the 
flow solutions. In order to achieve that, it will first be necess ary t o define a new physical variable, ip = p^'''^^'^^^vr'^ , closely following 
a perturbative procedure prescribed by_ Petterson et alj 1 198(lh an d 'Theuns & David' j 19921) for sph erically symmetric flows, and applied 
successfully later to thin disc flows I Ray boOsT Chaudhury et al. ^2006: Ray & Bhatta charjee 2007). It is quite obvious from the form of 
equations ([7} and J14t . that the stationary value of il) will be a global constant, t/jq, which can be closely identified with the matter flux rate, 
within a constant geometrical factor. A perturbation prescription of the form v(r, t) — vo{r) + v'{r, t) and p(r, t) — po(r) + p'(r, t), will 
give, on linearising in the primed quantities. 



7 + 1 



£_ 

Po 



+ 



Vo 



(42) 



with ^p' being a linearised time-dependent perturbation about the constant matter inflow rate, tpo- It is significant that the foregoing expression 
for Ip' is free of the explicit presence of a. Linearising in p' and v' about po and t^o, respectively, in both equations (O and ( I13t , and expressing 
p' and v' separately in terms of ip' only, will ultimately lead to a linearised equation for the perturbation as 
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= 0, 



(43) 



which is an expres s ion that is exactly the same as what can derived upon perturbing the stationary solutions of conserved axisymmetric 
inflows I Ray 20031 : Chaudhury et al. 20061) . Another aspect of equation ( 143 1 is that its form has no explicit dependence on the potential 
that is driving the flow. This is entirely to be expected, because the potential, being independent of time, will appear only in the stationary 
background flow. Arguments regarding stability will, therefore, be more dependent on the boundary conditions of the steady flow. As the 
form of the equation of motion for the linearised perturbation remains unchanged even for a flow in a fractal medium, and as the phy sical 
boundary conditions are also not altered in this case, the general conclusions reached earlier regarding non-fractal axisymmetric flows dRav 
Eo"o3) . can be extended here, and it can be safely claimed that under all reasonable boundary conditions, both the transonic and subsonic 
solutions will be stable. The only difference that will arise will be due to a fractal scaling of the amplitude of the perturbation. Restricting 
the perturbation to be like a high-frequency travelling wave, and carrying out a Wentzel-Kramers-Brillouin analysis on equation ( 143b . will 
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show that the time-averaged energy flux in the perturbation, T, goes as V'(7^ JiavlbooiEhaudhurv et alJbooel) . Since i/)o oc m, going by 



equation ( 141b . one can argue that JT ~ M '^^ , something that gives the expected fractal scaling for the amplitude of the perturbation. 



8 CONCLUDING REMARKS 

Two salient facts have emerged out of the entire theoretical exercise carried out here: fractal flows can be equivalently modelled as an 
effective continuum flow with a diminished local density, and that fractal properties of the flow have an adverse bearing on the question of 
multitransonicity. It is conceivable that there should be some observational imprints of these two aspects of fractal flows. The mass accretion 
rate in an accretion flow is intimately related to the local density field. It has already been argued here that the mass accretion rate undergoes 
fractal scaling, and there are observational indications of such scaling for the particular case of accretion processes onto a PMS star. 

In the context of multitransonicity, shocks are a very closely related feature. Many earlier works have discussed the transition of the 
flow from one transonic solution to the other via a shock (Chakrabarti 1990, 1996). Now fractal flows have been shown to be detrimental to 
multitransonicity. So, if shocks do or do not arise in a fractal flow, then that should cast some light on the actual mechanism behind shock 
formation in an accretion process. Again some observational signature of this should exist. 



The exac t manner in which matter infall takes place along the radius of an accretion disc, has been a subject of much study dLvnden-Befl 



I969I: |Pringlel[l98ll : iNaravan & Yilll994l:lFrank et al.ll2002h. The standa rd view is that viscosity (in some form) aids the outward transport of 



iBular momentum ('Shak ura^^unvae vn9730Balbus & Hawlevll998 ), and effects the formation of a Keplerian distribution of the accreting 



matter (jPringle 1981; Frank et al. 2002). What has been shown here is that for a conserved low angular momentum sub-Keplerian disc, a 



fractal structure in the accreting medium would drive the flow closer to the simple iBondil ( 119521) accretion limit, and would therefore assist 
in the infall process onto a black hole. So a combination of fractal behaviour and viscosity in the flow could possibly explain the observed 
infall rates. 
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